Astron. Nachr. / AN , No. , 1 - 1 () / DOI please set DOI! 



Interaction of radio jets with clouds in the ambient medium: Numerical 
simulations 
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Hydrodynamical simulations of jets interacting with clouds moving in the ambient medium of the host galaxy are pre- 
sented. Clouds with sizes of the order of the jet diameter and smaller, crossing the path of the jet with different speeds 
are considered. In the case of slow moving clouds the jet is stopped over the brief period of time taken by the cloud to 
cross the jet. The jet maintains its general morphology in the case of fast moving clouds. Erosion of the clouds leads to 
redistribution of cloud material to large distances. Such interaction may explain the large outflow velocities observed from 
pc to kpc scales in the compact radio sources. 
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1 Introduction 

The Gigahertz Peaked Spectrum ( GPS) and the Compact 
Steep Spectrum (CSS) sources (cf. lO'DealH^ . for a re- 
view) are young radio sources evolving in the Inter Stel- 
lar Medium (ISM) of the galaxies hosting them. The ra- 
dio morphology of these sources show evidence of inter- 
action of the radio jets with inhomogen eties in the ISM. 
( Jevakumar et al. I2QQ5I: Saikia et al ■Il995b . Several observa- 
tions of the line emitting gas in the optical and radio wave- 
lengths, show further evidence for such interactions. Obser- 
vations of the gas at large scales reveal that the line emis- 
sion from the narrow line regions (NLRs) and the cold neu- 
tral HI clouds show complex velocity profiles and outflows 



stand the relation between the outflow and the ionisation 
of the gas. In fact, several numerical studies have been at- 
tempted to study t he interaction of the jets with clouds in the 
ambient medium ( Choi, Wiita, & Ryu 2007;'Frasile et al.'2004'; | 
Higg i ns, O'Brien, & Dunlop 1999; Wang, Wiita, & Hooda 
|2QQQI: IXu & Ston jll995h . All these studies focused on the 
evolution of the individual clouds and the effect of the clouds! 
on the dynamics and morphology of the radio source. The 
clouds in these studies are large and massive resembling dif- 
fuse neutral clouds in the ISM of the Galaxy. A few simu- 
lations of jet propagation in a clumpy medium have been 
carried out where a few of them also considered parameters 
similar to that of the CSS/ GPS sources, but with a lower jet 



kinetic power (iSaxtiraet al. 2005; Steffen et al. 1997; Sut herland & Bickn elH 



dde Vries et al.ll999l:lEinonts et al.l2005l:lLabiand200"8[lLabia i>2m7iL| 

2005 ; 'Morganti. Tadhunter. & OosterloQl2005l: [o: Pea et al] In all of these simulations, the jet collides with a static 
2003). In addition, outflows have been observed in the broad cloud or a dense clumpy structure in the ambient medium, 
components. These observations reveal outflows of line emit-| However, a more realistic scenario is dense small clouds 
ting gas with velocities ranging from a few 100 km/s to such as those expected in the NLR and BLRs that move into 
more than 1000 km/s, from pc to kpc scales, both in the to the path of the jet/lobe. Such a scenario is more proba- 
ionised and the neutral components (cf. Holt, Tadhunter & ble for the young radio sources where the radio jet has not 
Morganti 2003, 2008; Morganti et al. 2005). Modelling of cleared the natal cocoon of material of the ISM yet. So, here 
lines of high ionisation state (Gu pta, Srianand & Saiki a 2005 ),| simulations of clouds moving into the jet are attempted, 
as well as several narrow lines (ILabiano et al.ll2005h . sug- 



gest that there is a contribution by shocks to the ionisation 
of the gas. All these observations suggest that interaction of 
jets with clouds in the ISM could play a role in the ionisa- 
tion and outflow of the gas. 

Numerical simulations have been successfully used to 



riety of ambient atmospheres (cf. Carvalho & O'Dea 2002; 


Hooda & Wiital 1 


9961: 1 Jevakumar et al.i.2005: Krause 2005; 


Wiita & Norman 


1992|). Such studies can be used to under- 



2 Jet cloud interaction 

The line emitting regions of the ISM are expected to be 
clumpy and turbulent. It can be approximated as a hot ten- 
uous inter-clump medium with spherical clumps embedded 
in it. There are two ranges of densities and scale sizes of 
the clumps in the ISM that show outflows: (a) the NLRs and 
dense large neutral clouds (n ~ 10^~^ cm~^) at kpc scales 
and, (b) BLRs with very high densities (n ^ 10^^ cm 
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but at sub-pc scales (cf . Ferlan d et al.ll992l:lLaor et al.l l997l;l 
Sulentic. Marziani. & Dultzin-Hacvanl2000l:lxu et al.l2007l) .| 
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Fig. 1 Snapshots of the logarithm of density (top) and the Lagrangian tracers of the cloud and jet material (bottom) of 
the run Fl are shown. The images from left to right correspond to 1.0, 1.04 and 1.12 problem time units. Cyan color (light 
gray) represents the cloud material and magenta (dark gray) represents the jet material. The radius of the jet and the size of 
the cloud are scaled to 100 pc. In this unit the total extent of the visible jet is about 3 kpc (left panel). 




Fig. 2 Same as Fig[Tl but for the run S 1 (slow cloud). The snapshot times are 1.0, 1.18 and 1 .65 problem time units from 
left to right. The radius of the jet is scaled to 100 pc and the size of the cloud is 100 pc. In this unit the total extent of the 
visible jet is about 3 kpc (left panel). 



To simulate a cloud with a typical density of BLR, a grid 
resolution of milli parsec is required. For the purpose of 
studying the overall dynamics involved in the interaction of 
the jet with a moving cloud, small dense clouds at two dif- 
ferent velocities are considered. The density of the clouds 
considered here is comparable to that of the NLRs. The 
sizes of the clouds are chosen such that at least 10 zones 
span the cloud diameter. Because of this constraint the clouds| 
are very massive. A smaller cloud spanning only 5 zones are 
also considered to reduce the mass of the cloud. 

2.1 Numerical setup 

Hydro-dynamical simulations in 3D were carried out using 
the ZEUSMP-V2.0 code, making use of the parallel fea- 



ture that has become available with this code ( Clarke|[r996l: 
iHaves et al.l l2006: St one & Norman 1992a,b). The jet prob- 
lem wasjnitiali^ed^ the previous jet launching code 
(cf. I Jeyakumar et al.ll2005D . For the simulations presented 
here a cluster of 20 processing cores were used. These sim- 
ulations were run on a Cartesian grid with 180 x 126 x 22 
active zones spanning 50 Rj in the jet propagation direc- 
tion, X, and ±15 Rj in the cloud propagation direction, y, 
and ±1.1 Rj in the z direction, where Rj is the radius of 
the jet. This choice allows us to consider a better resolution 
in two directions at the cost of a smaller physical extension 
in the third direction, although the resolution is the same. 
In the y direction, 80 uniform zones were used to span the 
inner 8 Rj region, while the rest of the zones were loga- 
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rithmically scaled. In the z direction all zones were uni- 
formly spaced. In the x-direction logarithmic scaling was 
used. This enables us to track the propagation of the cloud 
with a reasonable resolution. Outflow boundary conditions 
were used everywhere except at the jet inlet, where inflow 
boundary condition was used. 

The unit of length in this problem is the radius of the 
jet, Rj. The unit of time is given by the sound speed of the 
ambient medium, which was set to 1.0. The important pa- 
rameters of the jet are the intrinsic Mach number, Mj, and 
the density, t], with respect to the ambient medium density. 
The initially pressure matched jet is assumed to be conical 
with a half opening angle of Vl. The ambient density is given 
as, Tiair) = ^0/(1 + {r/a)'^Y, where a is the core radius, 
6 is the power law index of the density profile and r is the 
distance from the nucleus of the source. 

The cloud is characterised by the central density, ric, 
size, Rc, and a velocity with respect to the sound speed of 
the ambient medium, Vc. The density of the cloud is tapered 



as na{r) + {ric — na{r)) tanh((i?c — ^)/(0.3 Rc)), to avoid 
numerical artifacts (cf. Choi et al. 2007), where R is the dis- 
tance from the centre of the cloud. The temperature of the 
cloud is varied in such a way that the pressure inside follows 
the pressure of the undisturbed ambient medium, although 
the cloud is introduced inside the region affected by the 
bow shock of the jet, leading to a small but negligible pres- 
sure difference between the cloud and the disturbed ambient 
medium. In addition, Lagrangian tracer variables are used to 
track the jet and the cloud materials using the multi-species 
advection scheme available with ZEUSMP V2.0 (see Hayes 
et al. (2006), for more details). These variables are set as 1.0 
at the corresponding grid locations at the initial times. 

The numerical parameters used in these simulations are, 
Mj = 26, 7? = 0.001, a = 20 R^-, no = 1.0, S = 0.75 
and ft = 0.02. After the jet has propagated to about 30 Rj, 
the cloud is introduced at x=2 Rj, y=2 Rj and z=0. The 
numerical runs Fl and F2 (called fast cloud) with an Rc of 
1.0 and 0.5 Rj respectively are calculated with a velocity Vc 
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of 50. The numerical runs SI and S2 (called slow cloud) are 
calculated with the same parameters but with slVc of 5. 

These simulations are scaled as follows: for a jet radius 
of 100 pc, and temperature of the ambient medium of 2 x 
10^ K, the time unit is 0.43 Myr and the jet power is 2 x 
Iq4:5 QYg/s. For these parameters the velocities of the fast 
and slow clouds are 10^ and 10^ km/s and the core radius 
of the density profile of the ambient medium is 2 kpc. The 
central density of the cloud is 10^ cm~^. 

3 Results and Discussion 

Color images of the logarithm of the density as well as the 
Lagrangian tracers of the jet and the cloud are shown in 
Fig [T] to Fig m These figures show the morphology of the 
jet and the cloud at different times: (a) when the cloud is 
launched, (b) when the cloud is half-way through the jet 
and (c) when the cloud has just passed through the jet. Dur- 
ing the initial interaction of the cloud with the jet, erosion 
of the cloud material is seen (see the centre panels of the 
figures). The eroded cloud material is carried to large dis- 
tances. Apart from this, the cloud remains as a clump block- 
ing the path of the jet. At later stages, distortion of the cloud 
can be seen and the cloud material is distributed to a dis- 
tance of about half the length of the lobe. At later stages 
when most part of the cloud has passed through the jet, a 
clump of material can be identified as a cloud, but com- 
pressed in the direction of propagation of the jet. 

The time taken by the slow, small cloud to travel a dis- 
tance of a jet diameter is about 3 x 10^ yr, while the larger 
cloud takes somewhat longer time to cross through the jet. 
The jet is stopped completely in the slow cloud runs and the 
morphology is similar to that seen in other simulations of 
static clouds. In the fast cloud simulations, the cloud veloc- 
ity is 10 times faster, thus the interaction times are smaller. 
These simulations also show the general morphology seen 
in the slow cloud runs. However, the jet is not completely 
stopped in these runs and the radio source maintains its gen- 
eral morphology. Although the parameters of the jet and 
those of the clouds used here are somewhat different be- 
tween the simulations of interaction of the jets with static 
clouds (e.g. Higgins et al. 1999; Wang et al. 2000), these 
results can be compared with that of the static cloud simu- 
lations. The main difference is the erosion of the cloud ma- 
terial at the initial stages. Since the clouds are very massive, 
the clouds do not gain any velocity in the x direction. How- 
ever, the eroded material seen at large distances during the 
initial encounter accelerates to large distances and can be 
compared with the outflow. At later times (the right panels 
of the figures), back flowing cocoon also carries the mate- 
rial backwards. In such a scenario the emission lines can be 
very complex. 

In young radio sources where the radio jets are still ex- 
panding within the natal cocoon of a clumpy ISM, such in- 
teractions are likely to occur. Once the jet has emerged out 
of the ISM clearing the cocoon of material on its way, such 



incidences are less likely to occur. Obs ervations also indi- 
cate such a scenario (I nskip et al. I l2002h . These simulations 
also suggest that the distorted radio morphology and out- 
flow can be correlated. 
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